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If two elastic bodies with rough surfaces are first pressed against each other and then loaded tangentially, 
sliding will occur at the boundary of the contact area while the inner parts may still stick. With increasing 
tangential force, the sliding parts will expand while the sticking parts shrink and finally vanish. In this paper, 
we study the fractions of the contact area, tangential force and tangential stiffness, associated with the 
sticking portion of the contact area, as a function of the total applied tangential force up to the onset of full 
sliding. For the numerical analysis randomly rough, fractal surfaces are used, with the Hurst exponent H 
ranging from 0.1 to 0.9. Numerical simulations by boundary element method are compared with an 
analytical analysis in the framework of the Greenwood and Williamson (GW) model. In both cases, a 
universal linear dependency between the real contact area fraction in stick condition and the applied 
tangential force is found, regardless of the Hurst exponent of the rough surfaces. Regarding the dependence 
of the differential tangential stiffness on the tangential force, a linear relation is found in the GW case. For 
randomly rough surfaces, a nonlinear relation depending on H is derived. 

Surface roughness plays a widespread role in many physical phenomena including friction, wear, sealing, 
adhesion, electrical and thermal conductivity as well as in the acoustic properties of contact interfaces. The 
main basic contributions to the understanding of the role of surface roughness in tribological contacts are 
due to Bowden and Tabor\ Archard^ and Greenwood and Williamson (GW)l In the last decade, the contact of 
rough surfaces became once again a hot topic, and the focus moved to the analysis of fractal, self-affme surfaces 
showing roughness in a wide range of wave vectors. In particular, the real contact area A between rough surfaces 
was shown to be almost exactly proportional to the normal load F/'^ in a wider range of normal forces than 
predicted by the GW theory 

A(f.) = K^, (1) 

where Wz is the rms slope of the surface and k is a dimensionless factor. The reduced modulus of elasticity for two 
bodies in contact made of the same material is = E/{1 — v^) = 2G/(1 — v), with G being the modulus of shear 
and V the Poisson ratio. The normal contact stiffness of rough surfaces has been shown to be proportional to the 
normal load in the case of "nominally flat" surfaces (surfaces whose roughness power spectrum has a long 

wavelength cut-off or roll-off)^. For surfaces without a long wavelength cut-off, the contact stiffness was found 

1 

to be a power function of normal load^'^: kz oc F^^", where z is the coordinate axis perpendicular to the surface and 
H is the Hurst exponent. Tangential contact of rough surfaces was studied in^'^° and^\ In particular, in^ it was 
shown that the ratio of normal and tangential stiffnesses in a complete stick-contact remains the same as in the 
case of single rotationally symmetric contacts (Cattaneo and Mindlin ratio^^'^^). W\ this result was confirmed 
with surfaces of GW type and the tangential stiffness was found to rise linearly with the normal load in that case. 

In^°, it was shown that in the case of partial slip, also the maximum tangential "pre-slip" distance is related to 
the indentation depth in the same way as in a single rotationally symmetric contact. 

So far, important properties of rough contacts that are connected with the stick and slip regions of a tangential 
contact have not been investigated in detail. Yet the separate investigation of the stick and slip parts of the real 
contact area is relevant for many physical applications. For example, the differ ential tangential stiffness is only 
dependent on the current configuration of the stick region of contact. If bodies are subjected to tangential 
oscillation with a small amplitude (smaller than the maximum slip amplitude up to complete sliding), then wear 
will occur only in the slip-region (fretting). As shown by Ciavarella and Hills^^ for arbitrary two-dimensional 
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Figure 1 | Schematic drawing a single spherical aperity contact loaded tangentially. Due to the loading, the outer ring of the contact area experiences 
some slip. Displacements are exaggerated for visuaUsation purposes. 



contacts and in^^ for three-dimensional axis -symmetrical contacts, 
during fretting, the initial stick region always stays in the stick state 
(both at a constant normal force and a constant indentation), while 
in the initial sliding region the surfaces finally loose contact comple- 
tely. This leads to changes in stiffness as well as electrical and thermal 
conductivity, which in the final state will only depend on the initial 
configuration of the stick region. The area of the stick region is also of 
interest, for example for impurity-film resistance of contacts. The 
decomposition of the contact area into stick und slip parts plays an 
important role in further practical problems such as frictional energy 
dissipation^^'^^, geotechnical engineering^^, damping and stiffness of 
joints^^, and control engineering^^. 

In the present paper, we investigate the dependence of the sticking 
part of the contact area on the normal and tangential force, as well as 
its contribution to the total tangential force and to the total tangential 
stiffness. The investigation is based on the classical solutions for 
tangential contacts going back to the works by Cattaneo^^ and 
Mindlin^^, which later have been generalized to arbitrary rough con- 
tacts by Ciavarella^^'^^ and independently by Jager^^. They have 
proven that the solution of the tangential contact problem can be 
gained from the linear superposition of two normal solutions, one 
corresponding to the limiting frictional traction distribution and a 
negative correction due to a distribution equal to the coefficient of 
friction multiplied by the normal contact pressure distribution cor- 
responding to some lower value of the normal load. Numerical 
implementation of this technique was discussed in^^ and its extension 
to a more general loading case consisting in an oscillatory tangential 
force is available in^^. To date, although the mathematical theory and 
the numerical methods for the solution of this problem are estab- 
lished, the properties of the tangential contact solution for rough 
surfaces have been only marginally analysed. For smooth surfaces 
like spheres or cylinders, it has been shown in^^ that the tangential 
force component acting in the stick region Fx,stick and the real stick 
area Astick are both nonlinear functions of the total tangential force 
F^. We will show theoretically and numerically how these trends are 
modified in the case of rough surfaces. In particular, we will show 
that for the stick area a very universal linear dependence is valid, 
which practically does not depend on the character of roughness. 
When considering a tangential load that is monotonically rising, then 
the load has a biunique relation to the tangential displacement. In the 
numerical and analytical treatment of the problem, we will always 
impose the displacement and deduce the relevant quantities, and 
then express everything with respect to the tangential force at that 
point, in order to give an intuitive physical interpretation of the 
contact problem. 

Results 

Partial slip in multi-asperity modelling. For simplicity, consider a 
contact of an elastic half space having a flat surface with a rigid rough 
counterpart and assume that in each point of real contact the 



Coulomb's law of friction with a constant coefficient of friction pi 
is valid. We start our consideration with the basic version of the 
Greenwood and Williamson (GW) modeP and show later that its 
result is robust and does not depend on most of the assumptions on 
which the GW model is based. 

A nominally flat rough surface, consisting of a large number of 
independent "asperities", each having the same radius of curvature 
R, is brought into contact with an elastic half space. The final result 
does depend only very weakly on the exact form of the probability 
density function ^(z) of the height distribution. For simplicity, the 
probability distribution function is assumed to have an exponential 
form: 



0(2:) = Cexp(-2:/>l), 



(2) 



where is a characteristic length parameter of the height distribution 
("roughness"). 

First consider one single parabolic asperity subjected to a normal 
force Fzfi kept constant during the application of the tangential load, 
see fig. 1. Quantities named with '0' in the subscript refer to the case 
of a single asperity. 

The contact radius ao, the indentation depth d^, the radius of the 
stick region Cq and the tangential displacement are connected by 
the relation^^ 



Ux = fido - 



at 



(3) 



where, G* = 4G/(2 — v) is the reduced shear modulus and 
ao = \/Rdo applies. The maximum tangential displacement q up 
to the start of the full sliding can be found by setting Cq = 0 in (3) and 
is equal to 



--fido 



E* 



(4) 



At given Ux<uIq, full sliding takes place for all asperities for which 
0 < < d^m,o with 



dmin ,0 



fiE* 



(5) 



In contrast, if the indentation depth do of an asperity is larger than 
<^min,o» the asperity will be in the state of a partial slip. The normal 
force is given by the Hertzian equation^^ 



(6) 



The tangential force and the stick area for partial sliding are given 

by^^ 
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= HP..0 (l - I) = t^lE*R'/' [df - {do - d^^ ,of'] (7) 



and 

^stick,0 = TICq = nR(do — (imin ,o) • 

For full sliding we find 



and 



^stick,0 — 0. 



(8) 



(9) 



(10) 



Now we assume that the elastic half space is indented to the height 
z = h. The total stick area for a given tangential displacement is 
given by the integration of Eq. (8) over all the asperities that are high 
enough to have the indentation depth d^^^ necessary for achieving 
the partial sliding regime 



Astick= I nR{z-h-dmm)^(z)dz 

GO 

= nR^i^{i + h + d^i„)di 

0 

00 

= 7ti?C I ^ exp( - + /! + d„i„)/A)d^ 



(11) 



exp ( '-j^ I nRC ^ exp ^ 



Ao=Astick(Ux = 0 = dn,in) 



exp( -^1^0 

where Aq is the true contact area obtained from the normal contact. 
This corresponds also to the initial stick area for a vanishing tangen- 
tial displacement u^. 

The total tangential force is calculated by integration of (7) over all 
asperities in the state of partial sliding and of (9) over all asperities in 
complete sliding: 

GO 

F, = fi^E*R'/' I [{z-hf'-(z-h-d^^f']^(z)dz 

h + dmin 

+ fi^E*R'^^ I {z-hf^^(z)dz 

(12) 

GO 

= -^^E*R'/' I {z-h-d^^f'^(z)dz 

h + dmin 

GO 

+ fi^E*R'/^ ^{z-hf^^(z)dz 

h 

By substitution z — z and introduction of the notation 

GO 

F,,^,, = fiF, = fi^E*R'^' ^exp(-z/X){z-hf'dz (13) 



we can rewrite (12) in the form 

Fx = Fx,imx[l- exp(-^lniin/^)]. 

From Eq. (11) and Eq. (14) we finally get 



= 1- 



^stick _ ^slip 



(14) 



(15) 



Thus, for the classical GW modeP with an exponential distribution of 
asperity heights, we find linearity between the tangential force and 
the stick contact area. It can be shown easily that this result does not 
depend on the assumed paraboUc shape of the asperities. 
Assumption of any other asperity form or the consideration of an 
ensemble of fractal rough asperities leads to the exact same result. 
Note that this result does not depend on the length parameter X 
either. This fact implies a weak sensitivity of the solution with respect 
to the form of the probability distribution function. The only prop- 
erty needed for the derivation of Eq. (15) is a rapid (exponential) 
decay of the probability function. 

Furthermore, Eq. (15) can also be obtained for many more cases 
using the analogy of the tangential and normal contact problem. 
Indeed, it follows from^^ that in a tangential problem with partial 
slip, one can express the tangential stress Xx in the stick region as the 
sum of a full- sliding term and a correction term ql 



(16) 



where p is the pressure distribution. It shall be noted that q^/ fiis a 
solution for pressure to the normal contact problem. This corrective 
solution has the initial sticking area as its contact area and 



if 



<5* {x,y)dxdy 



(17) 



is the corresponding normal force. Integrating (16) over the whole 
contact area we find 



fx 



(18) 



When we are interested in quantities closely related to the stick 
region, it is sufficient to use the alternate normal contact solution 
at normal force = Q*. This theorem holds whether the stick region 
is a compact area or split into multiple fractions. It can be shown that 
the indentation depth in this corrective normal contact problem 
reads 



(19) 



The major quantities related to the region of stick are the area and the 
tangential stiffness of this region. Solving for the area of stick and 
assuming the validity of Archard's law (1) we can write 



Astick{Fx)=A{Q'') = K—^-K 



E*Vz fiE*Vz 



Astick{F,)=A{F,)(l-^^ 



(20) 



thus finding a very general formulation of eq. (15) for all systems, 
where Archard's law holds. These can be like those statistically 
described by GW, but also other nominally flat surfaces or fractal 
rough surfaces. 

For the normal contact stiffness as a function of the applied load, 
let us assume a general power-law dependency with exponent a, 
which is also expected from dimensional analysis considerations 
based on incomplete self-similarity^°: 
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Figure 2 | Numerical dependency of the normalized stick area as a 
function of the normalized tangential load. Bars denote a confidence 
interval of 95% with respect to the average result. For all values of the Hurst 
exponent H, the dependency is very close to the linear relation, Eq. (15). 



(21) 



For any (mostly isotropic) contact region, the tangential stiffness is 

(22) 

and we find for the tangential stiffness of the stick region 



k=^k 



(23) 



For fractal rough surfaces with Hurst exponent H and no cut-off or 
roll-off in the PSD"'': 



a = (H+l)" 
so we find in dimensionless variables 



= 1 



(24) 



(25) 



In the case of nominally flat surfaces, the PSD must be limited in the 
long-wavelength-part, e.g. by either a cut-off, or a roll-off (constant 
PSD). The latter case formally corresponds toH=— l,soa~l (see 
reference 8). The GW model also describes a particular type of nom- 
inally flat surface and for both cases we expect a linear dependency 

^-1-4- (26) 

It shall be noted that even though we can state general properties 
about the corrective stresses inside the stick zone, no general conclu- 
sion can be drawn regarding the stick force, because the pressure 
distribution inside the stick region is unknown. 

Nevertheless, we can do this kind of analysis for the GW model in 
case of an exponential distribution of asperity heights, where all 
quantities can be determined in closed form. Indeed, the integrals 
from Eq. (12) can also be separated in order to distinguish between 
the tangential stick force and the slip force. The slip force includes a 
component from asperities in partial slip and another from full- 




Figure 3 I Numerical dependency of the normalized differential 
tangential stiffness as a function of the normalized tangential load. Bars 
denote a confidence interval of 95% with respect to the average result. 
Numerical data support the dependency (25) shown in dotted lines. For H 
= 0.1, the greatest discrepancy is found. This is likely to be related to the 
fact that for small H, eq. (24) depends on the grid resolution, see [8]. 



slipping asperities. After some lengthy manipulation that we omit 
here for the sake of brevity, we find 

Fstick 



exp(- 



l)+erfc(^V^) 



(27) 



I^Fz 

where b = d^UX. 



= 1 — 2 exp( - 



/ 71 — erfc yVbj . 



(28) 



Numerical (BEM) results for fractal surfaces. Using the boundary 
element method (BEM) applied to numerically generated fractal 
rough surfaces, we find the stick area to follow (15) very closely, 
see Fig. 2. This result holds for all values of H that we investigated 
in the range from 0.1 to 0.9. For each value of H, 20 surfaces are 
generated and the bars of dispersion are plotted over the average 
response, to provide information about the scatter in the 




Figure 4 | Dependency of the normalized tangential load carried only by 
sticking regions as a function of the normalized overall tangential load. 

Numerical results are given for three values of the Hurst exponent H, with 
confidence interval of 95%. The solid green line shows the analytical 
dependency resulting from GW model, Eq. (27) together with Eq. (14). 
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Figure 5 | Evolution of the stick/slip contact domains obtained using BEM for an increasing tangential force. Black areas denote stick regions, grey is for 
sliding regions, white means no contact. With increasing tangential force (left to right), slip regions grow from regions supporting lower normal forces up 
to the whole contact area. 



numerical results. It is very interesting to see linearity especially for H 
~ 1 with a very small deviation from it. This result underlines once 
again that the distribution of surface height is irrelevant to the 
problem, because in the case of fractal surfaces it is far from being 
exponential (see reference 31 for a statistical analysis of the height 
field of surfaces generated using the random midpoint displacement 
method). Second, the elastic interactions neglected in the original 
GW theory used here and accounted for using BEM certainly play an 
important role on the traction distributions for fractal surfaces with 
large H, as also demonstrated for the normal load scenario in^^. Thus, 
numerical results vahdate the relation (15). 

The dependence of the normalized differential tangential stiffness 
on the normalized force is shown in Fig. 3. While for small Hurst 
exponents in the vicinity of H = 0, the dependence is close to the 
linear relation (26), for higher Hurst exponents, the GW Equation is 
not supported by numerical results, but approaches eq. (25). 

The fraction of the stick component in the generated tangential 
force is shown in Fig. 4. For small tangential forces, almost the whole 
contact domain is in the stick-state, so the stick force is almost 
identical to the overall tangential force. With increasing tangential 
load, more and more contact spots slip, causing the relative stick 
force to decrease. When the tangential force reaches its maximum 
value given by the normal force times the coefficient of friction, the 
stick areas vanish and so does the stick tangential force. Surprisingly, 
the numerical results for the surfaces with H = 0.9, which by defini- 
tion have strong asperity interactions, provide the closest trend to the 
GW prediction (27) plotted with solid line in the same figure. 

An example of evolution of the stick and slip contact domains is 
shown in Fig. 5. 

Discussion 

To conclude, we found a linear relation between the stick area Astkk 
and tangential force for nominally flat rough surfaces in the frame- 
work of the GW contact theory and for general tangential contacts 
where Archard's law holds. The results are quite universal and are 
completely supported by direct numerical simulations of contacts 
with self-affine rough surfaces with Hurst exponents in the interval 
from 0.1 to 0.9. We further investigated the differential tangential 
stiffness both theoretically and by means of direct BEM simulations. 
The theory predicts a linear dependency for nominally flat surfaces 
and H ~ 0. Indeed, numerical results show this linear dependency in 
the limit of small Hurst exponents and transition to a power-law 

behaviour k^/ k^^o = {l—Fx/ fiFz)'^ that is closely related to the nor- 
mal contact response. This result is obtained assuming the validity of 
the Ciavarella/Jager-principle and applying eq. (21), (24). 

Finally, we investigated the stick component of the tangential force 
and we found characteristic bell- shape curves for both the numerical 
samples and for the GW theory. The maximum value of these curves 
for fractal surfaces depends on the Hurst exponent H. 



Methods 

BEM calculations. The 3D Boundary element method was used to solve the problem 
of tangential contact with partial slip. As a first step, the normal contact must be 
solved using an iterative scheme. As a result, the apparent area of contact is 
partitioned into two sections. Discrete grid points that are considered to be in contact 
are all deflected to the same height and have positive pressure. All other points are not 
in contact, so they have zero pressure and a positive gap width. We used the CG- 
approach proposed by Polonsky and Keer^^ 

In the next step, a tangential displacement Ux is applied. This time, the true contact 
area from the first step is again divided into two parts. On the one hand there can be 
points that are in stick state, which means their local tangential displacement is equal 
to Ux and their tangential stress must be smaller that the local normal pressure times 
the coefficient of friction. On the other hand, points can be in sliding state, so their 
tangential stress is equal to the coefficient of friction times the normal pressure. A 
correct partitioning is found iteratively. In order to do so, we start from an arbitrary 
partitioning and we apply the tangential stress x = |ip in the slip zone. After that, the 
tangential stresses in the stick zone are solved demanding that the tangential dis- 
placement here shall be constant Ux. In case any grid point from the slip zone was 
deflected more than Ux or any point from the stick zone results with stresses greater 
than then these points are transferred to the other state. 

All discrete grid points influence each other via the elastic coupling, according to a 
discrete formulation of the Cerruti formula for a single force acting upon the elastic 
half-space. We used the Cattaneo-Mindlin analogy, assuming that the elastic coup- 
ling parallel to the surface can be expressed in the same way as normal deflections 
(Boussinesq), only modified by a coefficient of proportionality. We also performed 
BEM calculations considering the full anisotropic Cerruti solution and we get 
indistinguishable results. A comprehensive formulation of a partial- slip-algorithm 
can be found in^'*. 

We generated self-affine surfaces with Hurst exponents from 0.1 to 0.9 using the 
random midpoint algorithm. For these surfaces we solved the partial- slip problem for 
a range of tangential displacements Ux and recorded the resulting tangential force and 
size of the stick area. The simulations have been carried out for square contacts with 
free boundary conditions outside the contact region. 
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